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Abstract 

Theoretical predictions for B decay rates are rewritten in terms of the Upsilon 
meson mass instead of the b quark mass, using a modified perturbation expansion. 
The theoretical consistency is shown both at low and high orders. This method 
improves the behavior of the perturbation series for inclusive and exclusive 
decay rates, and the largest theoretical error in the predictions coming from the 
uncertainty in the quark mass is eliminated. Applications to the determination of 
CKM matrix elements, moments of inclusive decay distributions, and the B — > Xsj 
photon spectrum are discussed. 



1. Introduction 

Testing the Cabibbo-Kobayashi-Maskawa (CKM) 
description of quark mixing and CP violation is a 
large part of the high energy experimental program 
in the near future. The goal is to overconstrain 
the unitarity triangle by directly measuring its 
sides and (some) angles in several decay modes. 
If the value of sin 2/3, the CP asymmetry in 
B ipKs, is near the CDF central value 1), 
then searching for new physics will require precise 
measurements of the magnitudes and phases of the 
CKM matrix elements. Inclusive B decay rates can 
give information on \Vcb\, iKsl, iVi^l, and |Vtd|. 

The theoretical reliability of inclusive measure- 
ments can be competitive with the exclusive ones. 
For example, for the determination of \Vcb\ model 
dependence enters at the same order of Aq(^j-,/to^ ^ 
corrections from both the inclusive semileptonic 
B XcCD width and the B — > D*eD rate near zero 
recoil. It is then important to test the theoretical in- 
gredients of these analyses via other measurements. 

The main uncertainty in theoretical predictions 
for inclusive B decay rates arise from the poorly 
known quark masses which define the phase space, 
and the bad behavior of the series of perturbative 
corrections when it is written in terms of the pole 
mass. Only the product of these quantities, the 
decay widths, are well-defined physical quantities; 
while perturbative multi-loop calculations are most 
conveniently done in terms of the pole mass. Of 
course, one would like to eliminate any quark mass 
from the predictions in favor of physical observables. 
Here we present a new method of eliminating in 
terms of the T(IS') meson mass j2]. 



2. Upsilon Expansion 

Let us consider, for example, the inclusive B 
X^ev decay rate. At the scale /i = nib, 
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The variable e = 1 denotes the order in the modified 
expansion, and [3q = 11 — 2n//3. In comparison, the 
expansion of the T(IS') mass in terms of mf, has a 
different structure, 
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where t = hi\^/ (mbasCF)] and Cf = 4/3. In this 
expansion we assigned to each term one less power 
of e than the power of as- It is also convenient 
to choose the same renormalization scale /i. The 
prescription of counting a" in B decay rates as 
order e", and a" in tot as order e"~^ is the upsilon 
expansion. It combines different orders in the as 
perturbation series in Eqs. (|l|) and (^), but as it 
is sketched below, this is the consistent way of 
combining these expressions. 

At large orders in perturbation theory, the 
coefficient of a" in Eq. (|l]) has a part which 
grows as Cn\(3^~'^ . For large n, this divergence 
is cancelled by the a""*"^ term in Eq. (^), whose 
coefficient behaves as (l/as)(C/5)n!/3Q~^ [3~5]. 
The crucial l/ois factor arises because the coefficient 
of a^^^ in Eq. (0) contains a series of the form 
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(^£n-i _)_^n-2 ^_ _ _|_ which exponciitiates for large 
n to give exp(£) = ^/ {rribasCF), and corrects the 
mismatch of the power of as between the two series. 

The infrared sensitivity of Feynman diagrams 
can be studied by introducing a fictitious infrared 
cutoff A. The infrared sensitive terms are non- 
analytic in A^, such as (A^)"/^ or A^" In A^, and arise 
from the low momentum part of Feynman diagrams. 
Linear infrared sensitivity (terms of order -s/A^) 
are a signal of Aqcd effects, quadratic sensitivity 
(terms of order A^ In A^) are a signal of Aqqj-, effects, 
etc. From Refs. it follows that the linear 

infrared sensitivity cancels in the upsilon expansion 
to order (probably to all orders as well, but the 
demonstration of this appears highly non-trivial). 

Thus, the upsilon expansion is theoretically 
consistent both at large orders for the terms 
containing the highest possible power of /3o, and to 
order including non-Abelian contributions. 

The most important uncertainty in this ap- 
proach is the size of nonperturbative contributions 
to 77ix other than those which can be absorbed into 
TOft. If the mass of heavy quarkonia can be com- 
puted in an operator product expansion then this 
correction is of order ^qqi^/ {asTUbY by dimensional 
analysis. Quantitative estimates, however, vary in 
a large range, and it is preferable to constrain such 
effects from data. We use 100 MeV to indicate the 
corresponding uncertainty. Finally, if the nonper- 
turbative contribution to T mass. Ax, were known, 
it could be included by replacing tox by mx — Ax 
on the left hand side of Eq. (||). 

There are three surprising facts that are either 
accidental or indicate that the nonperturbative 
contributions may be small: 1) applications in 
terms of the T(2S') mass give consistent results 
2) the D Xev rate in terms of the tp mass works 
(un)reasonably well [||; 3) sum rule calculations for 
e"'"e~ bb find that the IS b quark mass (defined 
as half of the perturbative part of 'Tix(is)) is only 
20 MeV different from tox(is)/2 0. 

3. Application 

Substituting Eq. (||) into Eq. (|l|) and collecting 
terms of a given order in e gives 
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The complete order result calculated recently |^ 
is included. Keeping only the part proportional 
to Po, the coefficient of would be —0.035. The 



perturbation series, 1 — O.llSe — O.OSle^, is better 
behaved than the series in terms of the 6 quark pole 
mass, 1 — 0.17 e — O.lOe^, or the series expressed 
in terms of the MS mass, 1 + 0.30e + O.lSe^. The 
uncertainty in the decay rate using Eq. (^) is much 
smaller than that in Eq. (|l^), both because the 
perturbation series is better behaved, and because 
mx is better known (and better defined) than 
TOfc. The relation between \Vub\ and the total 
semileptonic B —> X^eD decay rate is 
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The first error is obtained by assigning an 
uncertainty in Eq. ^ equal to the value of the 
term and the second is from assuming a 100 MeV 
uncertainty in Eq. (^). The scale dependence 
of \Vub\ due to varying /i in the range nib/ 2 < 
fi < 2mb is less than 1%. The uncertainty in 
Ai makes a negligible contribution to the total 
error. Although B{B X^ev) cannot be measured 
without significant experimental cuts, for example, 
on the hadronic invariant mass, this method will 
reduce the uncertainties in such analyses as well. 

The B — > X^ev decay depends on both rrif, and 
•nic- It is convenient to express the decay rate in 
terms of mx and Ai instead of nib and TOc, using 
Eq. (I) and 



nib — nic = niB — niD - 
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where niB — {iniB* + r7iB)/4 ~ 5.313 GeV and 
rfiD = (3mD. + mo)/^ = 1.973 GeV. We then find 
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where the phase space factor has also been 
expanded in e, and the BLM subscript indicates 
that only the corrections proportional to /3o have 
been kept. For comparison, the perturbation series 
in this relation written in terms of the pole mass 



is 1 - 0.12e - 0.07BLMe^ 
proportional to Ai^2, Eq- 



r. Including the terms 
implies 
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I Vcfcl = (41.6 ± 0.8 ± 0.7 ± 0.5) x 10 
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where 77qed ^ 1.007 is the electromagnetic 
radiative correction. The uncertainties come from 
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assuming an error in Eq. (^) equal to the term, 
a 0.25 GeV^ error in Ai, and a 100 MeV error in 
Eq. (H), respectively. The second uncertainty can 
be reduced by determining Ai (see below). 

The theoretical uncertainty hardest to quantify 
in the predictions for inclusive B decays is the size of 
quark-hadron duality violation. This was neglected 
in Eqs. (|^) and (R). It is believed to be exponentially 
suppressed in the mi, — > oo limit, but its size 
is poorly known for the physical b quark mass. 
Studying the shapes of inclusive decay distributions 
is the best way to constrain this experimentally, 
since duality violation would probably show up 
in a comparison of different spectra. The shape 
of the lepton energy [11-14] and hadron invariant 
mass jl^, ^|lj] spectra in semileptonic B XcgD 
decay, and the photon spectrum in i? — > X^j [17- 
21] can also be used to determine the heavy quark 
effective theory (HQET) parameters A and Ai. The 
extent to which these determinations agree with 
one another will indicate at what level to trust 
predictions for inclusive rates. 

Last year CLEO measured the first two 
moments of the hadronic invariant mass-squared 



(sh) distribution, {sh 



and {{sh 



subject to the constraint Ee > 1.5 GeV [ p2[ 
Here ruD — {mn + 3m_D«)/4. Each of_ these 
measurements give an allowed band on the A — Ai 
plane. Their intersection gives (at order a^) [|2^ 



A = (0.33 ±0.08) GeV, 
Ai = -(0.13 ± 0.06) GeV^ 



(8) 



This agrees well with the analysis of the lepton 
energy spectrum in Ref. [ p^ , although the order 
Aqqj) /ml corrections not included in these analyses 
introduce large additional uncertainties |l4| |l6| . f 

In the upsilon expansion A is not a free param- 
eter, so we can determine Ai directly with smaller 
uncertainty. Considering the observable |l^ i?i = 
E,{dr/dE,)dE,/ {dT/dE,)dE,, a 

fit to the same data yields H 



Ai = (-0.27 ± 0.10 ± 0.04) GeV^ 



(9) 



This is in perfect agreement with the value of Ai 



implied by the CLEO result for {sh 



the 



upsilon expansion. The central value in Eq. (|9|) 
includes corrections of order a'^Po [^i but the 
result at tree level or at order as changes by less 
than 0.03 GeV^ The first error is dominated by 

t CLEO studied moments of the lepton spectrum pd] , but 
the band corresponding to (Ee) = (1.36 ± 0.01 ± 0.02) GeV 
on the CLEO plot cannot be reproduced using the formulae 
in Ref. [|l l| . So I consider only the result in Eq. (p|) from 
Ref. |22r (1 thank Iain Stewart for checking this calculation.) 
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Figure 1. Prediction for {1 — xb)\xb>i-s in the 
upsilon expansion at order (e^)BLM (thick solid curve) 
and e (thick dashed curve). The thin curves show the 
contribution of the O7 operator only. (From Ref. ]l^.) 



I /ml corrections [Q not included in Eq. (||). We 
varied the dimension-six matrix elements between 
±(0.5 GeV)"^, and combined their coefficients in 
quadrature. The second error is from assuming a 
100 MeV uncertainty in Eq. (||). 

Another way to test the upsilon expansion, or 
determine the nonperturbative contribution to mx , 
is from B Xgj- Possible contributions to the 
total rate from physics beyond the standard model 
are unlikely to significantly affect the shape of the 
spectrum. The upsilon expansion yields parameter 
free predictions for moments of this distribution. 
Experimentally one needs to make a lower cut on 
Ey, so it is most convenient to study 
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where xb — 2E^/mB- The parameter 5=1 — 
2E™"^/mB has to satisfy 6 > Aqcd/^ib, otherwise 
nonperturbative effects are not under control. 
Order a^Po corrections to the photon spectrum 
away from its endpoint were computed recently [p^ . 
Fig. 1 shows the prediction for {1 — xb)\xb>i-s 
as a function of S, both at order e and (e^)BLM, 
neglecting nonperturbative contributions to mx- 
A +100 MeV contribution would increase (1 — Xb) 
by 7%, so measuring (1 — a;^) with such accuracy 
will have important implications for the physics of 
quarkonia as well as for B physics. 

For E^ > 2.1 GeV Fig. 1 gives {I - xb) = 
0.111, whereas the central value from CLEO [^ 
is around 0.093. Interestingly, including the CLEO 
data point in the 1.9 GeV < E^ < 2.1 GeV bin, 

the experimental central value of {1 — xb) over the 
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Decay 
widths 


Expansions 

pole 1 

m|J and 


in terms of 
niT and e 


B -* Xcev 
B Xnev 
B X^TV 

B ^ XnTV 


l-0.12e-0.07e2 
l-0.17e-0.13e2 
l-0.10e-0.06e2 
l-0.16e 


l-0.10e-0.03e2 
l-0.12e-0.03e2 
1 - 0.07e - 0.02e2 
1 - 0.08e 


B — > X(.u(s + d) 
B Xcc(s + d) 
B ^uu{s + d) 
B X.uc{s + d) 


1 - 0.05e - 0.04e2 
l + 0.20e + 0.15e2 
1-O.lOe 
l + 0.09e 


l-0.03e-0.01e2 
l + 0.16e + 0.07e2 
1 - 0.05e 
l + O.Ue 



Table 1. Comparison of the perturbation series for 
inclusive decay rates using the conventional expansion 
and the upsilon expansion [^. The second order terms 
are the BLM parts only. 



region > 1.9 GeV is 0.117, whereas the upsilon 
expansion predicts 0.120. Ultimately, one would 
like to see whether prediction and data agree over 
some range of the cut i?™™. One can also evaluate 

(1 — xb) in terms of A and Ai without using the 
upsilon expansion. The CLEO data in the 
region > 2.1 GeV implies the central values 
Aa. ~ 390 MeV and A„2^„ ~ 270 MeV at order 



and a^Po, respectively 
not dominate at order a1 



The BLM terms may 
l2l||, so it is important to 



calculate the complete two-loop correction to the 
O7 contribution to (1 — xb). 

The upsilon expansion has been applied to form 
factor ratios in exclusive semileptonic B decays, 
as well as nonleptonic decays Q, where it also 
improves the perturbation series (see Table 1). 
However, the semileptonic B branching fraction 
or the average number of charm quarks in B 
decay agree with other predictions in the literature. 
Applications of similar ideas for e'^e" — > tt are 
reviewed by Teubner at this Conference |24[| . 

4. Conclusions 

• Using Tot and the upsilon expansion, i.e., 
assigning order to the order a" term in B 
decay rates and e"~^ to the a" term in the 
perturbative expression for my, is equivalent 
to using a short distance b quark mass. 

• It improves the behavior of perturbation series 
for inclusive B decays, and eliminates mf, 
altogether from the theoretical predictions in 
favor of mx in a simple and consistent manner. 

• It may lead to smaller nonperturbative effects 
(to the extent this is reflected in the behavior 
of perturbation series). 



• The biggest uncertainty is the nonperturbative 
contribution to tot unrelated to the quark 
mass. It will be possible to estimate / constrain 
this from data in the future. 
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